A metric approach for sound propagation in nematic liquid crystals 
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In the eikonal approach, we describe sound propagation near to topological defects of nematic 
liquid crystal as geodesies of a non-euclidian manifold endowed with an effective metric tensor. 
The relation between the acoustics of the medium and this geometrical description is given by 
Fermat's principle. We calculate the ray trajectories and propose a diffraction experiment to retrieve 
informations about the elastic constants. 
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I. INTRODUCTION 

Number of systems in condensed matter physics can be 
described by the same mathematical structures as those 
describing the gravitational field in General Relativity: 
this is the core of the analogue gravity (or effective ge- 
ometry) programme In General Relativity, gravita- 
tion is accounted by the distortions of spacetime, which 
is modeled by a pseudo-riemannian manifold (H of sig- 
nature (-,+,+,+). The properties of this spacetime are 
encompassed in the metric tensor g (or simply metric) 
and other quantities derived from it, such as the Riemann 
tensor or the Ricci scalar. Free-falling particles follow 
trajectories corresponding to the geodesies (curves of ex- 
tremal length) of the manifold. In particular, light paths 
correspond to the null geodesies of spacetime 



ds 2 = g^dx^dx 11 = 0, 



(1) 



where ds 2 is the square of elementary length along the 
geodesic and the g^ are the covariant components of g 
(the Einstein summation convention is assumed, where 
greek indices refer to spacetime coordinates and latin in- 
dices are related to spatial ones). 

As light trajectories can also be curved when crossing 
a refractive medium, Fermat's principle can be stated in 
terms of an effective geometry 0, 9 ■ For a given refrac- 
tive index n(r), where r is the position vector of the light 
wavefront, the light trajectories are interpreted as the 
geodesies of a non-Euclidean space with the line element 
dS 2 



dS 2 



\r)(dx 2 +dy 2 + dz 2 ), 



(2) 



Despite ^ represents only a three-dimensional non- 
Euclidean manifold, the corresponding geodesies are 
equal to the null ones in a static four-dimcnsionalpseudo- 
riemannian manifold with the time coordinate [4{ . Since 
we will deal with static refractive indexes, whenever con- 
sidering the time, the time coordinate appears in the line 
element 



ds' 



-c 2 dt 2 



dY,' 



(3) 



' Corresponding author: crms@fis.ufal.br 



where c is the speed of the light in vacuum and dT, 2 rep- 
resents the spatial line element obtained by the method 
that will be exposed in this paper. For the study that will 
be developed here, c is the velocity of the sound in the 
material of interest when it is homogeneous and isotropic. 

In this paper, we will apply this gravitational anal- 
ogyto describe the sound propagation in liquid crystals 
[5H7|, materials raising a considerable interest as they 
are both theoretically challenging (81-Hoj and of practi- 
cal interest (PVA screens, q-plates or cellphones [111]). 
They are made of anisotropic molecules (disc-like or rod- 
like) that present meta-states between the crystalline and 
liquid phases. One of these meta-states is the nematic 
phase: the molecules' axis are oriented on average along 
a specific direction called director (represented by the 
versor n) whereas their centers of mass are randomly lo- 
cated. Once in the nematic phase, the refractive index, 
the sound velocity and other macroscopic properties be- 
come anisotropic. Moreover, the nematic phase can ex- 
hibit birefringence, as the director plays the role of the 
optical axis [f| H EH • 

The gravitational analogy will be based on Fermat's 
principle, once it can also be used to determine sound 
trajectories [HI . The local speed of sound plays the role 
of the speed of light in General Relativity and sound 
propagates along the null geodesies of the acoustic met- 
ric. The interest of using this analogy (or, more properly, 
this geometrical approach) is that it helps understanding 



2 



from a different point of view several problems dealing 
with acoustics (for example, tomographic image recon- 
struction via ray tracing [14j . flight time in acoustical 
waveguides [IH, problems using variational method [Hj]) 



121- 



or nematoacustics 

We will consider the sound propagation around a punc- 
tual (called hedghog) and a linear (called disclination) 
topological defect of the nematic phase of a liquid crys- 
tal with rod-like molecules O]. Sound trajectories 
will be exhibited and the diffraction patterns will be de- 
termined in the perspective of retrieving the elastic con- 
stants of the liquid crystalline medium. 



N 2 sin 2 a 

Vp 

c77 



N 2 cos 2 a 

vp 



IK 



(5) 



Following Q , we define the refractive index vector N p = 
^ and denote by index j_ (resp. //) components of vec- 
tors that are orthogonal (resp. parallel) to director h. 
Thus, the components of N p are N p± = N p sina and 
Np,. = N p cos a. For acoustic waves in anisotropic crys- 
tals (which is the case of a liquid crystal in its nematic 
phase [2(], [22| ) , it is known that 



II. ANALOG MODEL OF ANISOTROPIC 
LIQUID CRYSTALS 



• k = 



N p = v. 



(0) 



Identifying (|6]) with ([5]), we found the components of iT g 



A. Acoustic aspects 

The specific structure of the liquid crystal in nematic 
phase allows two kinds of acoustic waves to exist [El.l6l.fl9l]. 
similarly for sound in solid crystals [2(| Hlj: the ordi- 
nary wave, that behaves as inside an isotropic medium, 
and the extraordinary one that depends on the angle be- 
tween the direction of the propagation and the director 
h. For each kind of wave, one usually defines two char- 
acteristic velocities: the phase velocity, which points in 
the direction of the wave vector k, and the group velocity 
that is oriented similarly to the acoustic Poynting vector 
S (in general, not parallel to k) and indicates the direc- 
tion of energy propagation. The phase velocity v p of the 
extraordinary wave, for a liquid crystal in the nematic 
phase with the director h orientated on the z-direction, 
can be measured by the pulse-superposition technique 
[lH, producing the expression 



N„ 



N 



v 2 p (a) 



(Cn + iC^-Cn) cos 2 a) 



(4) 



where p is the density, a is the angle between the wave 
vector k and n 5,6]. This relation expresses that a liquid 
crystal can be viewed as a simple solid with linear clastic 
constants in the x- and y- directions both equal to C\\ 
and in the z-dircction equal to C33 0, [l9| . Additionally, 
the smallest portion of liquid crystal where the director 
can be defined presents local cylindrical symmetry and 
therefore, we will consider Q to be locally valid for any 
configuration of the director n (even in the presence of 
defects as in sections IIII Al and MI B[) . 

The procedure to obtain N g (f3) relies firstly on the de- 
termination of v p (a) and phase refractive index N p (a). 

To determine the latter, one replaces (0| in N 2 = v 'i^ , 

where v is the velocity of sound in the isotropic phase of 
the liquid crystal, obtaining 
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J 9/l 



Vp 
C33 



(7) 



Introducing /3 as the angle between v'g and n, then alter- 
nate expressions for the group velocity components are 



v g± = v/NgSm/3 and v 
of ([7]) in © results in 



9/1 



v/Ng cos /3. Thus, the use 



»a ' N P 



W (N P± sin/3 



N p// cos 



0) 



N 2 



-^sin^+^cos 2 /? 



Finally, we obtain 



N 2 (P) 



^sin 2 /?- 



— cos 2 , 
033 



(8) 



In the next section, we will discuss about the connec- 
tion between the acoustic group index and the effective 
metric experienced by sound in the nematics. 



B. Geometric aspects 

We intent to study the paths followed by the extraor- 
dinary ray. This can be achieved by applying Format's 
principle |23j , for which the path followed by sound be- 
tween points A and B is the one that minimizes the in- 
tegral 



F = 



Ng(P)dI, 



(9) 



where I is the arc length of the path. As pointed out 
by Joets and Ribotta [2J| for light inside an anisotropic 
crystal, the use of Fermat's principle to determine sound 
paths is equivalent to calculating the null geodesies of a 
manifold with line element dE 2 = . gijdx l dx^ , where 
this null geodesic is the curve that minimizes the integral 



N 2 = 
a 1 p 



9 

vp 



Cn sin 2 a + C33 cos 2 a 



(10) 
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The mathematical resemblance of the two previous ideas 
in the application of © and (|10p suggests one to make 
the identification 



N 2 (/3)dl 2 =J29ijdx i dx j = dT? 



(11) 



where dl 2 is recognized as the Euclidean line element. 
Still according to Joets and Ribotta, the resulting space 
is a Finslerian one. However, due to the symmetry of 
the liquid crystal molecules and the spatial configuration 
of the director n in the applications where this metric 
approach will be applied, the Riemannian geometry is 
enough to describe the sound propagation [25|, H(| . 

Thus, substituting (fTTl) in we found the general 
expression of the line element experienced by the sound 
wave in cylindrical coordinates 



ds z 



c 2 dt 2 



+ N 2 ((3) (dr 2 +r 2 d6 2 + dz 2 ) 



(12) 



and the covariant components of the acoustic metric re- 
lated to the line element (TT2]) are 



V 



o 









r 2 N 2 ((3) 








N 2 (P) J 



(13) 



Form the standpoint of sound, an anisotropic medium 
can thus be mimicked by an effective gravitational field. 
As for light, the WKB (or short wavelength) approxima- 
tion enables to identify sound paths to the null geodesies 
of the effective metric g. In pseudo-Ricmannian geome- 
try 0, 0, , the curve that minimizes the line element 
ds 2 is obtained from the geodesic equation 



d 2 x fl „ dx^ dx a 

1_ r£ = 0. 

dr 2 va dr dr 



(14) 



with r a parameter along the geodesic and the com- 
ponents of the Riemannian connection 



dx v dx a dx^ 



(15) 



In the next section, we will apply these equations to 
study sound propagation near defects occurring in the 
nematic phase of liquid crystals. 





FIG. 1. Hedgehog and (k = 1, c = 0)-disclination defects. 
The arrows represent the director h. 



of the nematic phase: thepunctual defect called hedgehog, 
with director h = f written in spherical coordinates, and 
the linear defect called disclination with director h = p 
written in cylindrical coordinates (Fig. [TJ. 

To determine the acoustic metric associated to an 
hedgehog defect, the strategy is to express cos /3 and sin /3 
in terms of the laboratory frame coordinates, by the Eu- 
clidean line element dl 2 = dr 2 +r 2 (dd 2 + sin 2 6*d</> 2 ) , sub- 
stitute them in © and we obtain the g^s by ([TT|) . 

In Frenet-Serret frame |28j, the position vector along 
the sound path R = rf and the tangent vector T(l) are 
related by 



j<n\ ^ d(rr) dr „ ^ dr 
dl dl dl dl 

As T(l) has the same direction as v* g , then 

R » r, dr 

T ■ n = cos p = — = r. 

dl 

From the modified Euclidean line clement 



(16) 



(17) 



1 = r 2 + r 2 2 + sin 2 Qd> 2 



one identifies 



sin /3 = \ r 2 I 9 2 + sin 



(18) 



The replacement of (fTTf and (TT8"|) in © results in the line 
element 



III. APPLICATIONS 



A. Sound paths near a defect in nematics 



which gives the effective metric for sound in the vicinity 
of a hedgehog defect. The rescaling of the radial coordi- 
nate of this equation by r = . / jr—vr produces 



When a liquid crystal transits from the liquid (or 
isotropic) phase to nematic one, defects (points or lines) 
can arise spontaneously, causing a reorientation of the di- 
rector. We will focus on two different topological defects 



dY? = dr 2 + b 2 ? 2 (d9 2 + sin 2 



(19) 



where b 2 = ^ a . Therefore sound paths are the geodesies 
of (fB3)) , that is the spatial part of line element of a global 
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FIG. 2. Sound trajectories on the equatorial plane of a hedge- 
hog in a liquid crystal with b = 0.9. This result is similar to 
sound trajectories on z = C st plane near to a (k = 1, c = 0)- 
disclination lying on the z direction with b = 0.9. 



Here, we denote by = -^r, g^ are the components of 
the metric g (determinant g) given by eq. ([T^]). We have 
also used the convention of the repeated indexes for sum- 
mations. Solutions are the usual harmonic plane waves 
of the form $ = $(t, r, 9, <p) = e~ wt ^(r, 6, </>), where w is 
the wave pulsation. 

The spherical symmetry of the problem allow us to 
write i/j(r,0,^>) = 4>(r,d) = J2^Zq o-iRi(r) Pi (cos 9), where 
Pi(cos#) are the Legendre polynomials of order I and 
ai are worthless constants for the partial wave method. 
Applying ip(r,6) in $(t,r,6,cj)) and expanding (f2"2"j) . we 
obtain the radial equation 



Ri(r) 



w a_ 1(1 + 1) 



b 2 r 2 



Ri(r)=0, (23) 



monopole [29) . and they are depicted on Fig. © for 
b = 0.9. 

We now study the case of disclinations. These defects 
are expressed in terms of parameters k and c Q, which 
are related to director in cylindrical coordinates by 



ft = (cos (k<f> + c) , sin (k<f> + c) , 0) . 



(20) 



Repeating the same procedure as for hedgehog defects, 
we obtain a generalized effective metric, similar to one 
found in M 



dY>' 



P 

C33 

z. 

C33 



cos 2 a 



sin 2 a 



P 

Cu 



sin 2 a 



cos 2 a 



dr 2 



(21) 



sin a cos a 



rdrdcj), 



where a = k<j> + c. In the case of (k = 1, c = 0)- 
disclination, the geodesies are represented on Fig. [2] 

Equations (|T9|) and (f2T|) are examples of effective met- 
rics experienced by the sound close to topological defects 
in nematic liquid crystals. In the next section, they will 
serve as basis for calculations beyond the WKB approx- 
imation to determine scattering sections. 

B. Acoustical diffraction by defects in nematics 

The connection between geometric and wave optics is 
well-known: in the WKB approximation, light rays iden- 
tify with curves that are tangent at each point to the 
Poynting vector [22j ■ Similarly, since the direction of the 
energy velocity for sound waves defines the direction of 
the sound rays [2l| , we use the partial wave method [30[ 
to examine the scattering of sound waves by the previ- 
ous hedgehog defect and indicate its diffraction pattern 
(for disclinations, a similar analysis can be found in [3lj). 
Sound plane waves obey d'Alembert's wave equation for 
scalar fields: 

1 



'd u $) = 0. 



(22) 



where Ri(r) = d r Ri(r). Solutions to this equation are 
Bessel functions of the first kind Ri(r) = J„(n(r), where 



When b = 1, we recover the flat space and the last equa- 
tion becomes 1 + i. The phase shift 61(b) of the scattered 
wave is 



m = l(i+\-n(i) 



n 
2 



1 + 



1 + 



.(24) 



The angular distribution of the scattered sound is given 
by [3^ the differential scattering cross section a (6) 



a(0) = \f(9)\ 2 = 



1 OO 

— £(22 + l)(e M <-l)P ; (cos( 



(25) 



where f(6) is called scattering amplitude. The spherical 
symmetry of the scattering generates an annular diffrac- 
tion pattern, a ring of sound. The angular location of the 
maximum of sound is given by eq. (|25j) , as it is shown in 
Fig. ©. 

An analytical expression for the location of the light 
diffraction ring scattered by a global monopole defect in 
the real spacctime is found in [32j and a similar one is 
developed here. The idea is to expand the scattering 
amplitude about b 2 ~ 1 



/(0) = / (o) (0) + / (1) (0) 



(26) 



and to analyze the angular behavior of the first two 
terms. For the phase shift ([24]) . we make £ 



1+ i 



and a 2 = ^-r-, obtaining 



S,(b) 
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FIG. 3. Differential scattering cross section, eq. (|25[) . for a 
hedgehog defect with (b = 0.9, lj = 1) truncated at I — 600. 



We expand this last equation about a 2 ps (that is equiv- 
alent to expand about b 2 ~ 1), resulting in 



1 - 



2b( 



0(a 4 



(27) 



Substituting (|27l) in the scattering amplitude, we obtain 
the first two terms of the expansion, f^°\9) and f^'(9). 
They can be written in terms of the generating functions, 
h(8,a), of the Legendre polynomials 



h(9,a) 



E 

(=0 



e 7 ™^^ (cos 0) 



-\/2(cos7ra — cos( 



where a = 1 — 4. The zero-order term, f^ oS> (9), can be 
written considering the derivative of h(9, a) 



f { °\9) = 



sin 7TCV 



2\/2w (cos 7ra — cos 9) 3 / 2 ' 



(28) 



and the first-order term, f^'{9), can be written consid- 
ering h{6, a) 



2bu> ^/2(cos ira — cos ( 



(29) 



The equations (|28|l and ([2U)) diverge when 9 = 9q = net, 
that is the analytical expression for the angular loca- 
tion of the diffraction ring. For example, when b = 0.9, 
O = 7T (1 - 1/0.9) « 0.35 rad, as indicated by Fig. ©. 
For (k = l,c = 0) disclinations, the same result is ob- 
tained [3~0 ] and they are in agreement with the calcula- 
tions derived by Grandjean [33j ] - 



IV. CONCLUSION AND PERSPECTIVES 

We have developed an effective geometry approach 
to investigate properties of sound propagation inside 



anisotropic media such as nematic liquid crystals. In this 
framework, the influence of each defect is represented by 
an effective metric: metrics for hedgehogs are similar to 
those of global monopoles, whereas (k, c)-disclinations 
have more complex forms. Sound trajectories and diffrac- 
tion patterns are thus modified by the non-trivial met- 
ric. In particular, sound scattering by hedgehog defects 
is identical to light scattering by a global monopolc. 

It's interesting to note the relation between the met- 
ric approach shown here and Katanaev-Volovich's the- 
ory [34| . While Katanaev uses an affine transformation 
to deform and curve the elastic medium, generating the 
anisotropic properties, we start from the anisotropic ve- 
locity to derive the effective metric. 

We could rcinterpretate the obtained results by Cos- 
mology's view. Once sound waves can be understood as 
propagating perturbations in the effective metric gener- 
ated by the director field n(r), they are analog to gravi- 
tational waves propagating in a true gravitational metric 
generated by a cosmic defect . Thus, the annular diffrac- 
tion pattern due to defects in a liquid crystal can also 
be expected for gravitational waves coming on cosmic 
strings or global monopoles. In other words, the diffrac- 
tion pattern is a gravitational signature of the presence 
of cosmic defects, and its properties may help retriev- 
ing information about them. Furthermore, considering 
the weak effects of gravitational waves on possible detec- 
tors, we have presented an algebraic result for the angle 
of the maximum intensity of the scattered gravitational 
waves. Since the scattered wave on this angle has ampli- 
tude greater than incident wave's one, our results indi- 
cate the angular position where detectors of gravitational 
waves must be place to improve their power of detection. 
Another possible extension in this way is the possibil- 
ity to test a cosmological mechanis ms p redicted in the 
paper of T. Damour and A. Vilcnkin [35| : from our anal- 
ogy, production of sound waves should occur by cusps of 
disclinations in liquid crystals. If required by the referee, 
we can insert these comments in our article. 

In the language of cohomology [36| , the Volterra pro- 
cess for a cosmic string is obtained by a removing of a 
cylindrical solid angle, whereas the Volterra process of 
the corresponding antistring is obtained by an addition 
of a cylindrical solid angle. For example, in nematics, 
such pair will consist in a (k = 1, c = 0)-disclination, 
with b 2 < 1 in its effective metric, being the analogue 
for a cosmic string, and a (k = 1, c = 7r/2)-disclination, 
with b > 1 in its effective metric, being the analogue for 
the anti cosmic string. As it is known, the annihilation 
of disclinations in nematics, a velocity field is expected 
near the pairs, keeping the director fixed {371 - I39T ] . And, 
as in the analogue models based in moving fluid where 
the gradient of velocity field maps the gravitational field 
[4p| , we could previse the presence of a gravitational field 
in the vicinity of two annihilating cosmic strings. 

Still in cosmology, topological defects such as cosmic 
strings and global monopoles appear as a result of sym- 
metry breakings during phase transitions occurring while 
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the early universe cooled down. In the Kibble mechanism 
[4l| , the formation of these defects gives rise to emission 
of gravitational waves. In liquid crystals, defects appear 
as a result of SO(3) symmetry breaking of the isotropic- 
nematic transition phase during cooling. Formation of 
defects also obeys a Kibble mechanism 0, |42j ■ There- 
fore, by means of our analogy, we should be able (in 
principle) to characterize the properties of gravitational 
waves (such as polarization for example) from the study 



of sound waves propagating in nematic liquid crystals. 

Future developments of this work concern on the influ- 
ence of temperature on elastic constants, directors fluc- 
tuations (corresponding to fluctuations of the metric it- 
self (43J) and consequently on the sound trajectories and 
diffractions. Another interesting path should be the pos- 
sibility of simulating the exterior metric of spherical, 
uncharged, non-rotating mass, known as Schwarzschild 
metric [2|, |^ , for sound rays (as it was proposed for light 
in ft). 
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